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Equilibrium crystal surfaces, constrained to equilibrate by means of dissociative dimer deposition 
and evaporation, have anomalous global surface roughness. We generalize earlier results for one 
dimensional interfaces to two dimensions. The global surface width scales with surface size L 
as ~ log [I// (log L) 3] instead of the conventional form ~ logi. The surface roughening 
transition does not change in nature, but its location is subject to a large and slowly varying 
logarithmic finite-size-scaling shift. 
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I. INTRODUCTION 

Every now and then a new direction of research sheds 
an unexpected Ught on an older topic. In this instance, 
that older topic is equilibrium surface roughness, and 
a study of non-equilibrium driven surface growth trig- 
gered it. One dimensional (ID) surfaces lack equilibrium 
phase transitions, but while being driven their station- 
ary growing states can undergo roughening transitions. 
A recent example is directed-percolation-type roughen- 
ing 0]. While generalizing the latter to directed-Ising- 
type roughening (|J^, Noh, Park, and den Nijs discov- 
ered that the equilibrium point in their phase diagram 
had unusual properties. The surface width W, 



(1) 



a = ^ . This involves a representation of the equilibrium 
surface as a RW inside a N dimensional tube. In Sec. |^ 
we present our numerical data for a, Lx L lattice and gen- 
eralize the so-called healing time argument of Noh et al. 
from ID to 2D. Next, in Sec. 0, we map the 2D evenness- 
constrained equilibrium surface partition function onto 
that of a 2D surface with annealed random fields which 
couple to the surface heights, and use this to derive the 
relation Lfj-ee ~ L/{\ogL)i analytically. In Sec. VI we 



discuss the implications of the anomalous roughness to 
the roughening transition, and in Sec. VU we comment 



on possible ex perim ental realizations. We summarize our 
finding in Sec. 



VIII 



II. DISSOCIATIVE DIMER DYNAMICS 



X 



did not scale with the conventional (random walk) ex- 
ponent a = ^, but instead with what seemed to be a 
value close to a — ^. The origin of this is a global con- 
servation law, an evenness constraint on the occurrences 
of every surface height, due to the dissociative dimer de- 
position and evaporation dynamics. They related this to 
even- visiting random walks (RW) i.e., RW required 
to visit every site an even number of times. In follow-up 
studies 1^,^, it was shown that the value a = i is ex- 
act, using a mapping to Lifshitz tails |^] in the density of 
states of ID fermions in non-Hermitian random fields. 

In this paper we address the issue whether something 
similar happens in two dimensional (2D) surfaces. We 
show that the global surface roughness is again anoma- 
lous, i.e., that instead of the conventional logarithmic 
finite-size divergence, W'^ ^ ^ogL, the width diverges 
as ~ logLficc with an effective surface size L^cc ^ 
L/(logL)3. Moreover, the surface roughening transition 
temperature has a large and slowly converging logarith- 
mic finite-size-scaling correction. 

The outline of this paper is as follows. In Sec. || we re- 
view the most basic version of dissociative dimer dynam- 
ics leading to the global evenness constraint. In Sec. Ill 
we generalize the ID results of Noh et al. to semi-infinite 
N X L lattices and show that the width still scales with 



The dynamic process that gives rise to the global even- 
ness constraint has in its bare-bone version the follow- 
ing structure. Consider the so-called restricted solid- 
on-solid (RSOS) model on a ID or 2D lattice. Ev- 
ery lattice site contains an integer-valued height variable 
hr = 0, ±1, ±2, • • •, and nearest neighbors can differ only 
by Ah = 0, ±1. The latter means that only surface steps 
of height-one are allowed. The RSOS model has a long 
history in the theory of surface roughening transitions 
(see e.g., ref. g). 

Impose now the following dimer-type deposition and 
evaporation rule on the RSOS model. Particles can ar- 
rive on and leave the surface only in pairs (as dimers). 
Choose at random two nearest neighbor sites on the lat- 
tice. The surface heights of both sites are increased or 
decreased by one unit (with equal probabilities) but only 
when those two sites are at equal heights, and when the 
new configuration would not violate the RSOS rule (only 
single height steps). This simulates the deposition and 
the evaporation of horizontal dimers. The dimers are dis- 
sociative in the sense that the identity of a dimer is lost 
after it is deposited, and that a particle can arrive on and 
leave the surface with a different partner. 

This type of dynamics was studied recently in one di- 
mension PI and it was found that the equilibrium surface 
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width does not scale with a = i , as in conventional sur- 
face dynamics, but with the anomalous exponent o; = ^ 
instead. The value a = ^ is highly generic. Equilib- 
rium long range order cannot exist in ID, and therefore 
the probability to make an up or down step while walk- 
ing along an equilibrium surface is uncorrelated beyond 
a definite correlation length. This implies that the inter- 
face width scales like the dispersion of a random walker, 
i.e., that a = ^■ 

Dissociative dimer dynamics circumvents this generic 
picture by imposing a non-local constraint on the Gibbs 
distribution. Particles are deposited at each height level 
in pairs, and therefore the RW must visit each level an 
even number of times. This is a non-local constraint, due 
to the dissociative character of the dynamics. For exam- 
ple, on a locally flat surface segment, two dimers are 
able to land next to each other, but next, the two par- 
ticles in the middle are allowed to switch partners and 
evaporate together leaving two disconnected monomers 
behind. Repeated processes like this leave no local trace 
of the dimer nature of the dynamics, but globally ev- 
ery height level needs to be occupied an even number 
of times. In contrast, non-dissociative dimer dynamics 
obeys the constraint locally and therefore yields the con- 
ventional scaling exponent a = ^, just like monomer dy- 
namics. 



III. SEMI-INFINITE LATTICES 

We start with a generalization of the one dimensional 
1 X L lattice results of Noh et al. to TV x L lattices with 
N — 2, 3, • • •. On the one hand, having several chan- 
nels instead of one weakens the evenness constraint and 
maybe sufficiently to alter the a = exponent. On the 
other hand, critical exponents change typically only with 
dimension, and these N channel lattices are still ID sys- 
tems at large length scales. 

Our simulation results are shown in Fig. |. We plot 
the effective values of a as functions of the system size L 
ioT N — 2 and 3, evaluated as 



log [WiL)/Wi^L)] 
log 2 



(2) 



The roughness exponent remains close to a = ^ for both 
N = 2 and N = 3. The convergence is not impressive 
but comparable to that in the original one channel case. 
These large finite-size-scaling corrections are explained 
at the end of this section. 

The numerical results for a suggest we seek an ana- 
lytical generalization of the even-visiting random walk 
representation of the system and from that show that 
the surface roughness exponent does indeed not change. 
Consider a walk along the equilibrium 2 x L surface from 
one edge to the other (from x — to x — L). The x coor- 
dinate plays the role of time and going up or down along 
the surface is the RW aspect. The walk does not obey 



the evenness constraint along one specific channel, y = I 
or y = 2, because dimers can be deposited in two orien- 
tations, along and perpendicular to the channel. How- 
ever, if you walk twice, one time in each channel, then 
the constraint is still obeyed for the two combined walks. 
The evenness constraint is weakened, but obviously still 
present. 

On a formal level, the connection between the ran- 
dom walk and surface roughness follows from the well 
known transfer matrix formalism for evaluating the par- 
tition function of the surface. Let's first consider a ID 
lattice with monomer-type dynamics (no evenness con- 
straint). The labels n and uq in the partition function 
Z{n\no)x represent fixed boundary conditions of the sur- 
face at both ends, h — uq ed site 0, and h = n at site x. 
Z{n\no)x obeys the recursion relation 

Z{n\no)x+i = Z{n\no)x + Z{n - l\no)x + Z{n + l|no)a;, 

(3) 

which is a discrete version of the diffusion equation and 
can be reinterpreted as the time evolution of a ID ran- 
dom walker with Z{n\no)x the unnormalized probability 
to find the walker at position n at time t = x starting 
from position no at time t = 0. 

Let's generalize this to the two channel lattice, n now 
becomes a two-component vector, n = (711,712). The 
transfer-matrix-type recursion relation 



Z(ni,n2|no)a;+i 



E 



Z(ni - 7ni,n2 - m2|no)3;, (4) 



could be easily solved numerically because the transfer 
matrix is still finite. However, for the sake of the evenness 
constraint, we will pursue that Eq. (|^) resembles again a 
diffusion equation and again can be reinterpreted as the 
time evolution of a RW; but this requires some care and 
the introduction of an absorbing wall. 

This is a walk on a 2D lattice (711, 712). The dimension 
of the walker's space is equal to the number of channels, 
not to the dimension of the surface. h{l^x) represents 
the ni coordinate and /i(2, x) the 712 coordinate of the 
walker at time t — x. For example, when the walker is 
located at site (4, 5) at time t = 10, the surface height at 
site (1, 10) is equal to ft, = 4, and at site (2, 10) equal to 
h = 5. 

The integer-valued summation labels 77ii and m2 in 
Eq. ^ are subject to the RSOS condition that nearest 
neighbor surface columns can differ only by Ah = 0, ±1. 
This translates into two restrictions: Firstly, at every 
moment in time, the two coordinates of the walker can 
differ only by 7*2 — tii = 0, ±1. So the walk is restricted 
in the (ni, n2)-plane to a strip centered along the diago- 
nal. Til = 7i2, as illustrated in Fig. ^. In other words, the 
walk remains quasi one dimensional after all. Secondly, 
due to the same RSOS rule, each positional component 
ni and 712 of the random walk can change only by or 
±1 during each time step, t = x ^ x + 1. It is easy 
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to see that the walker can jump to only 6 new positions 
from any diagonal site, ni = n2, and to only 5 from every 
ofF-diagonal one, ni ^ 77.2 (in addition to not moving at 
all). 

In the equilibrium Gibbs distribution for the RSOS 
monomer-type surface dynamics, all surface configura- 
tions have equal probabilities. This is reflected in Eq. (Q) 
by the fact that all transition probabilities are equal. 
However, this recursion relation does not represent yet 
a proper Master equation for a random walk because it 
does not conserve probability. A Master equation needs 
to be of the generic form: 

Z(n|no)t+i - ^("'' n)^(n'|no)t, (5) 

and the transition probabilities, w(n',n), must conserve 
probability from any state n', i.e., 

Ywin',n) = l. (6) 

n 

The one channel transfer matrix equation of motion, 
Eq. (||) is of this form; except that we have to divide 
all transition probabilities by a common factor, equal to 
3. The two channel transfer matrix equation of motion, 
Eq. (I) does not satisfy Eq. (§). 

Partition functions scale in 'time' t = x exponentially 

as 

Z(n|no). ~ Ag, (7) 

with Ao the largest eigenvalue of their transfer matrices. 
In Eq. (H), Aq is simply equal to 3, i.e., the same com- 
mon factor that normalizes the transition probabilities. 
Such a simple common factor does not exist for Eq. (^) 
because the numbers of sites to which the walker can hop 
are different at diagonal and off-diagonal sites. This lack 
of conservation of probability can still be incorporated 
in the random walk representation by interpreting the 
boundaries of the strips as so-called absorbing walls. In 
the absence of the walls, the common factor would have 
been equal to 3^. The remaining exponentially decaying 
factor (Ao/3^)^ represents the absorption by the walls. 
Moreover, Ag-type factors drop out of the calculation of 
thermodynamic averages, which in the random walk rep- 
resentation is equivalent to performing all averages with 
respect to the so-called surviving ensemble only. 

In a free 2D random walk, the fluctuations in the two 
components decouple. They do so also in our random 
walk with absorbing walls on the strip, because we aver- 
age over the surviving ensemble only. This implies that in 
the direction parallel to the strip, the root-mean-square 
displacement scales with the conventional power as 

^\\it)-t-2, (8) 

at time scales much larger than a characteristic correla- 
tion time arising from the range of hopping distances of 



the RW during each time step. The root-mean-square 
displacement in the perpendicular direction, £,±{t) re- 
mains finite due to the presence of walls. Thus we re- 
cover (although in a somewhat contrived manner) the 
well known result that for a semi-infinite 2 x L lattice 
with L ^ 00, the equilibrium surface width still scales as 
in ID, with Vr= [(^j|(L)+ei(i))/2]5 - Li 

Let's return now to the reason for constructing this 
random walk interpretation, i.e., to deal with the even- 
ness constraint. For the one channel lattice, the core step 
towards the analytic derivation of a — ^ was the intro- 
duction of Ising variables to keep track of where the RW 
has been before. An Ising spin S'„ = ±1 was associated 
with each site, and flipped for each time interval the RW 
resided on that site. The evenness constraint was en- 
forced by requiring that at time t = x all spins point in 
the same direction as at time t = 0. We need only to 
establish a generalized equivalent formulation to demon- 
strate that the surface roughness exponent retains its ID 
value a = ^ on the two channel lattice. 

Introduce an Ising spin S{n) = ±1 to each diagonal 
site {n,n), and associate it with the plaquette consisting 
of sites {n — l,n), {n,n -f 1), {n 4- l,rt), {n,n — 1), and 
(n, n) on the strip (shown as the shaded area for n = 2 
in Fig. 1^) . This spin is flipped for each unit time interval 
the walker resides on (n — 1, n), (n,n + I), (n + 1, n), or 
(n, n — 1). S{n) does not flip when the walker resides 
on the diagonal site (n, n) because then the number of 
occurrences of height n increases by an even number, 
h{l,x) = h{2,x) = n. Note that two Ising spins flip 
when the walker occupies an off-diagonal site, i.e., S{n) 
and S{n + 1) for the site (n, n + 1). 

In the one channel case, the Ising spins act as gate 
keepers along the chain, and the RW cannot pass site n 
without fiipping its spin S{n). Similarly, in the two chan- 
nel case, the RW moves on the quasi ID strip and again 
a line of Ising spins act as gate keepers. Topologically 
and at large length scales, the situations are equivalent 
and therefore a must be the same as in ID. 

One detail seems different. On the chain, the walker is 
unable to pass site n without flipping its spin S{n), while 
on the strip, it can pass through plaquette n without 
flipping S{n). That might seem an essential difference, 
but it is not. A relaxed model where a spin is flipped 
stochastically when the walker passes was shown by Noh 
et al. in the context of the purely ID surface to belong 
to the same universality class as the deterministic spin 
flip model. This followed most clearly from the transfer- 
matrix equivalence to Lifshitz tails in the density of states 
of ID fermions in non-Hermitian random flelds. The ran- 
dom walker is the fermion and the Ising spins generate 
the random potential (see ref. |^ for details). 

The generalization to N x L surfaces is straightfor- 
ward. The RW moves in a dimensional tube (i.e., 
still remains effectively ID) and the evenness constraint 
still requires only a single line of Ising spins positioned 
along the body-diagonal of the tube. Each flips when 
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the walker comes within range; S{n) flips every time any 
coordinate of the RW position is equal to n. We con- 
clude from the above evenness-constrained RW construc- 
tion that the roughness exponent must indeed retain its 
ID value a = ^ for all finite N. 

Let's return to our numerical results. The shapes of 
the curves in Fig. ^ represent so-called finite-size-scaling 
corrections. These corrections are larger and more com- 
plex than for the purely ID one channel lattice This 
can be understood qualitatively as follows. For = 1, 
the evenness constraint is ineffective at small L such that 
aes (see Eq. (^)) decreases monotonically from a value 
near the free unconstrained RW exponent a = ^ at small 
L towards q = i at L — > oo. (A characteristic crossover 
system size L could be constructed in terms of the ra- 
tio between L and length scale L{^oo{L) over which the 
surface fluctuations are unimpeded by the evenness con- 
straint; Lfroo(i) is defined in the healing time argument 
of Sec. 0.) 

In contrast, the N — 2 and iV = 3 curves start off 
from a value near zero, overshoot a = i, and then bend 
backwards. At small L < N, the surface behaves as if 
it is 2D and as if the evenness constraint is absent, i.e., 
W ^ (logL)2. This explains why the curves start off 
at a ~ 0. The surface starts to behave one dimensional 
beyond L ~ iV, but initially remains still unaware of 
the evenness constraint, such that the a curves in Fig. |l| 
overshoot (towards a ~ i) and only back over at larger 
L where the evenness constraint kicks in. 



6(i)~(iiogi:)5. 



(11) 



This might seem at odds with the canonical Vl power 
implied by the RW interpretation; but only until one 
remembers that the walk is truly random only at time 
scales L much larger than a characteristic correlation 
time arising from the maximum hopping distance of the 
walker during each time step, and realizes that the lat- 
ter (projected along the tube's body-diagonal) is pro- 
portional to ^/N, and thus diverges simultaneously with 
L. For the sake of curiosity, we checked and confirmed 
Eq. (|ll|) numerically. Fig. || shows the average root- 
mean-square displacement, £_{L) ~ [(C(|(i) + £,'^{L)) / L]^ 
of the RW as a function oi L — N . 

Before we present our numerical results for the 
evenness-constrained surface width in 2D, we like to 
present an educated guess of what the behavior might 
be by generalizing the so-called healing length argument 
of Noh et at from ID to 2D. This argument was the 
least rigorous of their analytical results, but explained the 
anomalous exponent a = i at a simple intuitive level. In 
ID and also for all our semi-infinite N x L surfaces of 
the previous section, the argument runs (in a somewhat 
modified form) as follows. 

Consider the semi-infinite N x L surface. There must 
exist a crossover length scale Lfrco, larger than the sur- 
face width N but much smaller than the total channel 
length L, within which the evenness constraint can be 
ignored and the surface fluctuates freely such that the 
surface width increases as 



IV. ROUGHNESS ON L x L LATTICES 

We now turn our attention to dissociative dimer dy- 
namics on a truly 2D L x L lattice. For monomer-type 
dynamics, it is well known that the height-height corre- 
lation function diverges logarithmically. 



5(r) = {{hr 



-ro ^ro) 
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ttK, 



G 



log? 



and the surface width scales as 



1 
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logL, 



(9) 
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with Kc the so-called effective Gaussian coupling con- 
stant which in the RSOS model at T — s- oo (our numerical 
simulations) takes the specific value Kq — 0.9 The 
goal is to find out how the evenness constraint affects this 
on a global scale. 

In the random walk representation of the previous sec- 
tion, we need to take the limit where the dimension of 
the tube (N) and the random walk time (L) diverge si- 
multaneously, as = L ^ oo. This diminishes the 
usefulness of the RW representation. For example, ac- 
cording to Eq. (p^, the root-mean-square displacement 
of the RW along the body-diagonal of the tube must scale 
(for monomer dynamics) as 



I^(Lft.ec) - L 



free : 



(12) 



just like in ID monomer dynamics. The basic premise of 
the healing length argument is that at larger length scales 
the surface must be preoccupied satisfying the evenness 
constraint such that the surface width does not increase 
any further. Then, we need only to find out how Lf^cc 
scales with L. 

Divide the surface strip < 2: < L into segments of 
length Lfree. Choose one of these segments and mark 
all surface heights that are visited an odd number of 
times in that segment. It is known that these defect sur- 
face heights are uniformly distributed within the range 
\h\ < 0{W) such that their number is also proportional 
to . Suppose that the surface compensates for these 
defects one-by-one along the rest of the strip (ignoring 
that the additional surface fluctuations introduce new 
defects as well as repair simultaneously other old ones). 
The average root-mean-square distance the surface must 
fluctuate vertically to reach that specific defect height is 
again proportional to W, and that requires a typical hor- 
izontal length interval Ax ~ .Z^froc- So every segment of 
length Lfrcc leads to only one repair on average. Since 
there need to be W of them, it follows that 



Lb, 



(13) 
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and, using Eq. dH), that 



V. SURFACE ROUGHNESS IN RANDOM MEDIA 



Lr^WW\ and VF-L^, 



(14) 



such that for all x L lattices the roughness exponent 
is equal to a — ^. 

For the generalization of this argument to L x L lat- 
tices, we divide the 2D lattice in blocks of size LfrccXLfice- 
Within each block the surface fluctuates as if the global 
evenness constraint does not exist, and the surface width 
scales as ~ {271X0)"^ log Lfrcc- Assume that again, 
like in ID, the oddly-visited surface heights are uniformly 
distributed and therefore that their number is propor- 
tional to W. The same repairing scheme as above pre- 
dicts again one repair on average per block, i.e., that 



L 



-'free 



and therefore that 

In the limit of large L, we can invert Eq. (09) to 



W\L) 



2ttK( 



L 



(logL): 



(15) 
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(17) 



The crucial features are the logarithms and the j power. 
The former originate from the logarithmic divergence of 
unconstrained surface roughness in 2D and the latter 
from the scaling of the number of blocks, (L/Lfrcc)^, 
with dimension D. 

Compared to the conventional logarithmic divergence, 
Eq. (|l(]|), the lattice size L is replaced by a logarithmically 
modified effective length, Lfroe ~ L/(\ogL)i . Logarith- 
mic effects are notoriously difficult to confirm numeri- 
cally, except when you expect them. 

In Fig. ^, we present our numerical simulation re- 
sults. The equilibrium surface width for monomer 
and dimer dynamics are plotted on semi-log scales as 
functions of the system size L for L = 8, 16, 32, and 
64 (squares and circles, respectively). The slope of the 
monomer curve confirms the value Kq ~ 0.9 of Eq. (10). 



The dimer roughness line is also quite linear, but with 
a reduced slope, which means that its surface roughness 
could be fitted (in this range of length scales) by a simple 
logarithm, but with an enhanced effective Kq- However, 
when we plot versus L/W'^ as sugg ested by Eqs. (|l|) 
and (17), the slope of the (again straight) line is the same 
as in the monomer case, i.e., we regain the correct value 
of Kq. Therefore we conclude that for dissociative dimer 
dynamics, Eqs. ( p^ ) and Eq. ( p^ ) are indeed correct. In 
the next section, we put this result on a more rigorous 
analytic footing. 



Solid-on-solid models come in several variations. The 
RSOS model we used above, where nearest neighbor 
columns can differ only by Ah = 0,±1, is most con- 
venient for numerical simulations. The so-called discrete 
Gaussian model P,[Tot, where the A/i-restriction is re- 
laxed and replaced by a Gaussian interaction energy 



"^[^^1 - fcS^ = E 



(18) 
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between nearest neighbor columns, and h^' (the mean- 
ing of the ( ) brackets in the summation), is more ap- 
pealing in analytical discussions. In this section we use 
the latter. 

The global evenness constraint can be incorporated in 
the partition function in the following manner. Define 
integer- valued variables v^, one for every surface height. 
Notice that they are not associated with any specific lat- 
tice site. Each vt is equal to the number of times surface 
height level h appears in the specific configuration. The 
partition function, subject to the global evenness con- 
straint, can be written as 



CXp {-A[hr]) 



(19) 



with fugacity z = —1. Every configuration in which any 
of the Vh^s is odd receives zero weight. 

Next, we introduce randomness degrees of freedom 
c;i = 0, 1 associated with those same surface height levels 
and write the partition function as 



Z = ^ cxp 



{c,h} 



'hVh 



log 2) - A[hr 



(20) 



with z = cxp(— /i). Eq. ( |2C| ) can be interpreted as the 
partition function of an equilibrium surface in the pres- 
ence of annealed random external fields {ch} that sup- 
press the occurrences of specific surface heights. The 
noise is global in the sense that the random variable Ch 
affects surface height level h equally at every position r 
along the surface. For every occurrence of height level h, 
the Boltzmann weight is multiplied by a factor z when 
Cfi = 1 (and not if Ch = 0). The fugacity parameter 
z allows us to interpolate between monomer and dimer 
dynamics because at z = 1 we retrieve the conventional 
discrete Gaussian partition function. 

We are following in this closely the analogous formula- 
tion for the ID surface by Noh et al. . They found that 
in ID z = acts as a special point, a stable fixed point 
in the sense of renormalization transformations with all 
— 1 < z < 1 as its basin of attraction. We will show that 
this remains true in 2D. 

The limit z = is special. The fugacity Boltzmann 
factors become delta functions, z^''"'' S{chVh), such 
that the surface is strictly prohibited to pass through all 
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height levels for which Ch = 1. Each of them acts as an 
impenetrable barrier and the surface is restricted in its 
vertical height fluctuations to wander between two such 
randomly placed walls. The partition function factorizes 
and reduces to 



(21) 



where ^(2^) is the probability of finding two neighbor- 
ing Cfi = I walls at a distance Ah — 2^ apart. Walls 
are placed at random with probability ^, such that 
P(2^) = 2^24 2^ in Eq. (|2^) is the partition function of 
a surface restricted to fluctuate between two such hard 
walls: 



E 

{hr} 



n^(^-i'^ri) 



exp{-A[hr]). 



(22) 



The two following modifications cannot affect the large 
length scale behavior. Firstly, below the equilibrium 
roughening transition the integer-valued height vari- 
ables /ip = 0,±1,±2, • • •, can be approximated by real- 
valued continuous variables with the coupling constant 
K renormalized to its Gaussian value Kq. This is well 
known from the theory of conventional surface roughen- 
ing Secondly, we can replace the hard impenetra- 
ble walls by soft ones and approximate by 



(23) 
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where the trace runs over ordinary configurations of 
the Gaussian model but restricted to the subset where 
the root-mean-square surface width in each configura- 
tion is equal to ^. The evaluation of is elemen- 
tary. Compared to the Gaussian partition function, 
Zq — } 6xp(— A[ft,r]), it gains an exponential factor. 
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Zoexp [-aoL^ exp(-47rii'GC^)] , 



(24) 



where ao = jq. Finally, the surface width of the annealed 
ensemble average follows from Eqs. ( pl| ) and (^ using 
the method of steepest descent as 



e exp [-2^ log 2 - a„L^ exp(-47r JCgC')] 
^5 exp [-2^ log 2 - aoL^ expi-inKae)] 
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(logi)^ 
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This is the same scaling form as proposed in the previous 
section from the heahng argument, Eq. (^. 

The above derivation remains equally valid for small 
values of z around z — 0. There, the surface can tunnel 
through the walls, but with exponentially small prob- 
abilities (proportional to the length of the intersection 
contours) actually resembling more closely the soft walls 
in the above derivation. At the positive z side, it is in- 
tuitively reasonable that z = I (the monomer dynamics 
point) is the horizon and limiting point for this type of 
scaling. At the negative z side, the horizon extends to 
and apparently includes the point 2 = — 1. The numerical 
results and also the healing argument of the previous sec- 
tion suggest this. This is also intuitively consistent. Neg- 
ative values of z imply an imaginary chemical potential 
and therefore a mix of positive and negative Boltzmann 
weights that cancel out against each other; this enhances 
the exponentially decaying nature of the |z|'"'''^'' weight 
factors, and preserves this decay even at 2; = — 1. 



VI. ROUGHENING TRANSITIONS 



The surface is presumed to be above the roughening 
transition temperature Tr in most of the above discus- 
sion. For example, although Eqs. (|l^) and (|^ apply to 
all T, it is correct to approximate the discrete Gaussian 
height variables by continuous Gaussian ones only for 
T > Tn. Everywhere in the rough phase the global sur- 
face roughness has a logarithmic correction, Eq. (|T^), to 
its conventional simple logarithm. The natural question 
arises whether and/or how this influences the nature of 
the surface roughening phase transition. In this section 
we will show that the evenness constraint has no effect on 
the properties of the roughening transition in the thermo- 
dynamic limit, but gives rise to strong finite-size-scaling 
corrections, including an apparent shift in the transition 
temperature. 

Equilibrium surface roughness has been a topic of ac- 
tive research over several decades. The roughening is an 
example of a Kosterlitz-Thouless (KT) transition [Q. 
Exp erimental realizations include Helium crystal surfaces 
|]l5| , metal surfaces [|6| , and organic crystals |0 . 

The anomalous roughness exists only at a global scale. 
Within the crossover length scale Lfieo ^ L/{\ogL)i (de- 
fined in Sec. [V), the surface fluctuates freely in disregard 
of the global evenness constraint; i.e., inside the bulk the 
equilibrium surface remains indistinguishable from that 
in monomer dynamics, and since phase transitions are 
ruled by the bulk, the roughening transition must remain 
in the KT universality class. 

The 'evenness boundary effect' is, however, very 
strong. We will approach this from the perspective of the 
original argument by KT where they estimated the tran- 
sition temperature by the free energy of a single vortex in, 
e.g., superfluid films ll4|. The discreteness of the height 
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variables, hr = 0,±1,±2,---, in the discrete Gaussian 
solid-on-solid model plays the same role as those vortices. 
This might be somewhat surprising, but follows mathe- 
matically from a so-called duality transformation . In 
the rough phase, the heights can be treated as continu- 
ous variables and this represents the 'vortex- free phase'. 
Consider the so-called sine-Gordon correlation function 
in the Gaussian model, in which the /ir's are continuous. 



(exp [27ri(/i; 



r+ro 



exp 

-2-r/Kg 



1 



§(2^)'3(r) 



(26) 



(see also Eq. ^). The logarithm of this is the free en- 
ergy of placing two topological objects, exp(±27rift,), at 
a distance r apart. They favor integer values of h. That 
free energy scales logarithmically with distance r, and 
a single one costs therefore an amount of free energy 
/ = (tt/Kg) logi or / = (Tr/i^G— 2) logL if we allow it to 
be at any position in 2D space. The temperature Kq = ^ 
where the latter is equal to zero is the famous KT esti- 
mate for the transition temperature. (Recall that the 
coupling constant Kg is measured in units oi ksT such 
that Kg ~ T~^.) This is how this works for monomer 
dynamics. 

For dimer dynamics the two sine-Gordon 'charges' in- 
teract logarithmically, just like in Eq. (|2^) but only at 
distances smaller than Lfrco- Closer to the lattice size 
the interaction levels off such that a single excitation now 
costs 



/ ~ — log Lfrco - 2 log L. 



This yields, by setting / = and using Eq. (|l7|), 



Kg 



bo + log(logL)' 



logi 



(27) 



(28) 



with bo a constant. This is a significant finite-size- 
scaling correction, both in magnitude and in the slow- 
ness with which it scales with L. The transition appears 
to happen at a higher temperature. For example, at 
L = 1000 (larger than most experimental surface het- 
erogeneity lengths [|l6l), the shift is of order 



ATr log(logL) 



logL 



7%. 



(29) 



VII. EXPERIMENTAL REALIZATIONS 

The scope of this paper is foremost theoretical, i.e., 
to show how a seemingly benign topological conserva- 
tion law in the dynamics strongly affects the equilib- 
rium state. Experimental realizations of the evenness 



constraint should be possible to establish, but this re- 
quires careful considerations and collaborations with ex- 
perimentalists. For example, molecular-type bonded 
molecules like N2 look promising, but unfortunately, in 
solid N2 the molecules do not dissociate. As a starting 
point we include here a section with a brief discussion of 
the most essential features. 

Simple models like ours are meant and perfectly suited 
to discover fundamental scaling laws associated with sur- 
face growth, but they are rather simple minded when it 
comes to making direct contact with actual experiments. 
Those detailed theories of crystal growth have reached a 
high level of sophistication in recent years . The even- 
ness effect represents a topological feature which will be 
preserved with increased realism. After identifying a suit- 
able experimental system, the evenness constraint needs 
to be embedded into the appropriate more detailed the- 
oretical description. However, at this stage, it suffices to 
focus on general aspects, in particular those that possibly 
upset the conservation law. 

The condition in our model that the X2 molecules land 
only horizontally is not a serious constraint, because the 
molecule must dissociate as well. The latter requires that 
the binding energy between the two atoms be weaker or 
at best be of the same order of magnitude as the bond- 
ing energies inside the solid. In that case, both atoms 
strongly prefer to be close to the surface, and a vertically 
landed molecule will quickly decay. 

A more serious issue is surface diffusion, which is omit- 
ted in our model. The evenness constraint is preserved 
as long as a surface wanders around on the same terrace, 
but is lost when surface atoms jump across steps to lower 
or higher levels. Diffusion across steps is reduced by the 
presence of so-called Schwoebel barriers . If those bar- 
riers are strong enough, the jumps will occur infrequently 
enough that the evenness constraint remains satisfied at 
length scales large compared to the heterogeneity length 
scale of the experimental surface. The latter is the length 
scale at which defects like impurities pin and limit the 
surface dynamics; this length rarely exceeds lOOOA. 

Moreover, recall that the anomalous surface roughness 
is stable with respect to varying the fugacity parameter 
z in Eq. (^9|), i.e., it remains present when the evenness 
constraint is not strictly obeyed but only statistically. 
Therefore it might well be that diffusion across steps pre- 
serves the anomalous roughness beyond the above esti- 
mate. This issue needs further theoretical study. 

Finally, the search for experimental realizations does 
not need to be limited to dimer- type dynamics. The 
anomalous global surface roughness exists and is the 
same for X„-type dissociative dynamics with any n > 1. 
This was found to be true earlier in ID jsj, and follows 
also in 2D from generalizing our analytic arguments. 
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VIII. SUMMARY 

In this paper we studied how equilibrium surface 
roughness in two dimensional surfaces is affected by a 
global constraint that every surface height be present an 
even number of times in every configuration. We pre- 
sented numerical and analytic evidences. 

In semi-infinite N x L surfaces, the evenness constraint 
is weakened compared to the 1 x L surface, but not 
enough to change the global width. It still scales anoma- 
lously a,s W ^ just like in purely ID surfaces 

In truly 2D LxL surfaces, the global surface roughness 
is also anomalous, but the effect is weaker. The surface 
width scales logarithmically as 



W^L)^ log Lt,, 



(30) 



similar to conventional monomer dynamics, but with a 
logarithmically corrected effective surface size 



^free 



(logL)^ 



(31) 



The surface roughening transition does not change in na- 
ture, but its location is subject to a large and slowly 
varying logarithmic finite-size-scaling shift, Eq. ([29|). 

These results remain valid when the evenness con- 
straint is only statistically obeyed. Moreover, the an- 
alytic derivation in Sec. M involves a mapping onto a 
surface model with annealed randomly placed barriers, 
placed horizontally to the surface, that inhibit the verti- 
cal surface fluctuations. The surface roughness of those 
surfaces behaves the same as the evenness-constrained 
ones. 

Finally, 2D represents a critical dimension for ordinary 
equilibrium surface roughness. This follows trivially by 
evaluating the surface roughness in the Gaussian approx- 
imation, W ^ L", with a = {2 — D)/2 (logarithmic in 
D = 2). The surface width does not diverge in dimen- 
sions D > 2 and the surface remains asymptotically flat. 
Still, the healing length and random field arguments can 

2D 

be generalized to higher dimensions: L^cc L^+o and 
W ~ (ifree)" ~ . So for example, in £> = 3, 

monomer dynamics yields W ^ L~2 and dissociative 
dimer dynamics W ^ L~5. The evenness constraint al- 
ways flattens the surface at a global scale. 



ACKNOWLEDGMENTS 

This research is supported by the National Science 
Foundation under grant DMR-9985806, and by the Brain 
Korea 21 Project. 



[1] U. Alon, M.R. Evans, H. Hinrichsen, and D. Mukamel, 

Phys. Rev. Lett. 76, 2746 (1996). 
[2] J.D. Noh, Hyunggyu Park, and M. den Nijs, Phys. Rev. 

Lett. 84, 3891(2000). 
[3] H. Hinrichsen and G. Odor, Phys. Rev. Lett. 82, 1205 

(1999). 

[4] G.M. Cicuta and M. Contedini, ^nath. PR/9903063 

(1999); G.M. Cicuta, M. Contedini, and L. Molinari, J. 

Stat. Phys. 98, 685 (2000). 
[5] J.D. Noh, Hyunggyu Park, Doochul Kim, and M. den 

Nijs, Phys. Rev. E., in press (2001). 
[6] M. Bauer, D. Bernard, and J.M. Luck, J. Phys. A 34, 

2659 (2001). 

[7] R. Friedberg and J.M. Luttinger, Phys. Rev. B 12, 4460 
(1975). 

[8] M. den Nijs, J. Phys. A 18, L549 (1985). 

[9] S.T. Chui and J.D. Weeks, Phys. Rev. B 14, 4978 (1976). 

[10] J.V. Jose, L.P. Kadanofl, S. Kirkpatrick, and D.R. Nel- 
son, Phys. Rev. B 16, 1217 (1977). 

[11] H.J.F. Knops, Phys. Rev. Lett. 39, 766 (1977). 

[12] For reviews on the theory of surface roughness, see e.g., 
J.D. Weeks, in Ordering m Strongly Fluctuating Con- 
densed Matter Systems, edited by T. Riste (Plenum, New 
York, 1980), p. 293; H. van Beijeren and I. Nolden, 
in Structures and Dynamics of Surfaces, edited by W. 
Schommers and P. von Blanckenhagen (Springer, Berlin, 
1987); M. Wortis, in Chemistry and Physics of Solid 
Surfaces VII, edited by R. Vanselow and R.F. Howe 
(Springer, Berlin, 1988). 

[13] M. den Nijs, in The Chemical Physics of Solid Surfaces, 
edited by D. King and D.P. Woodruff (Elsevier, Amster- 
dam, 1994), Vol. 7, Chap. 4. 

[14] J.M. Kosterlitz and D.J. Thouless, J. Phys. C 6, 1181 
(1973); J.M. Kosterlitz, J. Phys. C 7, 1046 (1974). 

[15] RE. Wolf, F. Gallet, S. Balibar, and E. RoUey, J. 
Physique 46, 198 (1985); Y. Carmi, S.G. Lipson, and 
E. Pohurak, Phys. Rev. B 36, 1894 (1987). 

[16] For a set of recent reviews on metal surface properties, 
see e.g.. The Chemical Physics of Solid Surfaces, edited 
by D. King and D.P. Woodruff (Elsevier, Amsterdam, 
1994). 

[17] See e.g., R.F. P. Grimbergen, H. Meekes, P. Bennema, 
H.J.F. Knops, and M. den Nijs, Phys. Rev. B 58, 5258 
(1998). 

[18] Examples of reviews from the vast literature on epitax- 
ial surface growth; A.-L. Barabasi and H.E. Stanley, 
Fractal Concepts in Surface Growth (Cambridge Univer- 
sity Press, Cambridge, 1995); T. Halpin-Healy and Y.-C. 
Zhang, Phys. Rep. 254, 215 (1995); A. Pimpinelh and 
J. Villain, Physics of Crystal Growth (Cambridge Uni- 
versity Press, Cambridge, 1998); J. A. Venables, G. D. 
Spiller, and M. Hanbucken, Rep. Prog. Phys. 47, 399 
(1984); H. Brune, Surf. Sci. Rep. 31, 121 (1998). 

[19] G. Ehrlich and E.G. Hudda, J. Chem. Phys. 44, 1039 
(1966); R.L. Schwoebel and E.J. Shipsey, J. Appl. Phys. 
37, 3682 (1966). 



8 



0.45 



0.35 



0.3 



0.25 







0.035 
T -1 



0.07 



FIG. 1. The effective values of a in Eq. (|2|) v.s. the inverse 
of L are shown for 1 x L surfaces (diamonds), 2 x L surfaces 
(circles), and 3 x L surfaces (squares) with L = 16, 32, 64, 
128, and 256. 
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FIG. 2. The 2 X L (2 channel) equilibrium surface repre- 
sentation as a random walk on a strip with absorbing walls. 
The surface coordinate along the channel represents time in 
the random walk. The surface heights in the first and second 
channel represent the n\ and n2 coordinates of the walker. 
The shaded area is an example of a plaquette with Ising spin 
5(2) = ±1. The heavy solid lines are the absorbing walls. See 
the text for more details. 




FIG. 3. Semi-log plot of the mean-square displacement 
^^(L) of the random walk in the tube described in Sec. ^ as 
a function of time L for 2 < L < 40 with the total random 
walk time L equal to the dimension of the tubular random 
walk space iV. The sohd line fits the form ^{L) = 0.272 log L 
as predicted by Eq. (|ll|). 
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FIG. 4. Semi-log plots of equilibrium surface widths W 
on Lx L lattices for: monomer dynamics (squares) and disso- 
ciative dimer dynamics (circles) both as functions of L. The 
slopes of the straight-line fits yield the values of Kg (defined 
in Eq. (|l|)): Kg = 0.916 (monomer) and 0.988 (dimer). 
The diamonds show the same dimer dynamics data plotted 
as v.s. L/W^. That line has nearly the same slope as 
the squares one (monomer dynamics), i.e.. Kg ~ 0.914 in 
agreement with Eq. (p^. 
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